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Abstract. We infer the gas kinematics, diagnostics and ionic radial profiles, distance and central star parameters, 
nebular photo-ionization model, spatial structure and evolutionary phase of the Planetary Nebula NGC 6741 
by means of long-slit ESO NTT+EMMI high-resolution spectra at nine position angles, reduced and analysed 
according to the tomographic and 3-D methodologies developed at the Astronomical Observatory of Padua (Italy). 
NGC 6741 (distance~2.0 kpc, age~1400 yr, ionized mass Mi on —0.06 Mq) is a dense (electron density up to 12 000 
cm -3 ), high-excitation, almost-prolate ellipsoid (0.036 pc x 0.020 pc x 0.018 pc, major, intermediate and minor 
semi-axes, respectively), surrounded by a sharp low-excitation skin (the ionization front), and embedded into a 
spherical (radius~0.080 pc), almost-neutral, high-density (n(HI) ~7xl0 3 atoms cm -3 ) halo containing a large 
fraction of the nebular mass (Mh a i >0.20 Mq). The kinematics, physical conditions and ionic structure indicate 
that NGC 6741 is in a deep recombination phase, started about 200 years ago, and caused by the quick luminosity 
drop of the massive (M» =0.66-0.68 M Q ), hot (logT*~5.23) and faint (log L,/L Q ~2.75) post-AGB star, which 
has exhausted the hydrogen-shell nuclear burning and is moving along the white dwarf cooling sequence. The 
general expansion law of the ionized gas in NGC 6741, V cxp (km s -1 ) = 13xR", fails in the innermost, highest- 
excitation layers, which move slower than expected. The observed deceleration is ascribable to the luminosity 
drop of the central star (the decreasing pressure of the hot-bubble no more balances the pressure of the ionized 
gas), and appears in striking contrast to recent reports inferring that acceleration is a common property of 
the Planetary Nebulae innermost layers. A detailed comparative analysis proves that the latter (i. e. "U" -shaped 
expansion velocity field) is a spurious, incorrect result due to a combination of: (a) simplistic assumption (spherical 
shell hypothesis for the nebula), (b) unfit reduction method (emission profiles integrated along the slit), and (c) 
inappropriate diagnostic choice (A4686 A of He II, i. e. a thirteen fine-structure components recombination line). 
Some general implications on the shaping mechanisms of Planetary Nebulae are discussed. 
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1. Introduction 

The sentence by the late Professor Lawrence H. Aller 
(1994): "A nebula is a three-dimensional structure for 
which we obtain a two-dimensional projection" fully syn- 
thesizes the many, so far unsolved, observational lim- 

* Based on observations made with: ESO Telescopes at 
the La Silla Observatories (program ID 65.1-0524), and the 
NASA/ESA Hubble Space Telescope, obtained from the data 
archive at the Space Telescope Institute. Observing programs: 
GO 7501 and GO 8773 (P.I. Arsen Hajian). STScI is operated 
by the association of Universities for Research in Astronomy, 
Inc. under the NASA contract NAS 5-26555. We extensively 
apply the photo-ionization code CLOUDY, developed at the 
Institute of Astronomy of the Cambridge University (Ferland 
et al. 1998). 

Send offprint requests to: F. Sabbadin, sabbadin@pd.astro.it 



itations and interpretation problems connected to the 
Planetary Nebula (PN) research, leading to: (a) rough 
spatio-kinematical reconstruction, (b) unrealistic assump- 
tions for the gas parameters (in particular, electron tem- 
perature (T e ) and electron density (N e ) constant all across 
the nebula), and (c) proliferation of kinematical, physical 
and evolutional models, frequently based on the mere neb- 
ular morphology (i. e. appearance). 

At last, the tomographic and 3-D analyses developed 
at the Astronomical Observatory of Padua (Sabbadin et 
al. 2004 and references therein) overcame the stumbling 
block of nebular de-projection, rebuilding the ionic spatial 
structure, and allowing us a direct comparison of each real, 
true PN with the current theoretical evolutionary models 
(Schonberner et al. 1997, Steffen et al. 1998, Marigo et 
al. 2001), the detailed hydro-dynamical simulations (Icke 
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Fig. 1. Apparent morphology of NGC 6741 in the 
medium-excitation line A5007 A of [O III] (upper panel), 
and in the low-excitation line A6584 A of [N II] (lower 
panel). North is up and East to the left. HST-WFPC2 
frames; programme GO 8773, P.I. Arsen Hajian. 

et al. 1992, Frank 1994, Mellema 1997, Perinotto ct al. 
2004a), and the updated photo-ionization codes (Ferland 
et al. 1998, Ercolano ct al. 2003). 

Though the observational starting point is common - 
i. e. long-slit spectra -, the "philosophy" of tomography 
is just opposite of the conventional method. In fact, the 
latter compacts the spectrum along the slit (in practice, 
it restricts the nebula to a point), and gives mean, inte- 
grated results (line flux, expansion velocity, T c , N c , ion- 
ization etc.). Vice versa, tomography is based on a pixel- 
to-pixel analysis of both flux and velocity, and furnishes 
the bi-dimensional structure (in different ions) of the ra- 
dial slice of nebula intercepted by the spectrograph slit. 
Later on, a 3-D rendering procedure combines all tomo- 
graphic slices and provides the true spatial distribution 




Fig. 2. Apparent I([N II])/I([0 III]) distribution over 
NGC 6741 (same WFPC2 images as Fig. 1). The sharp, 
inhomogeneous layer at low-excitation ([N II]) framing the 
main nebula is punched, along and close to the apparent 
major axis, by a series of radial [O III] rays? yets? pene- 
trating into the faint, roundish halo. 

of the kinematics, physical conditions (T e and N e ), and 
ionic and chemical abundances at unprecedented accuracy. 
Tomography needs spectra at high "relative" spatial (SS) 
and spectral (RR) resolutions (SS=r/Ar, r=apparent ra- 
dius, Ar=seeing; RR=I4x P /AV, AV=instrument spectral 
resolution). It is based on the simple consideration that 
the position, depth and density of each elementary vol- 
ume within an extended, regularly expanding nebula can 
be, in principle, derived from the radial velocity, FWHM 
and flux, respectively, of the corresponding emission. 

So far we have studied NGC 40 (Sabbadin et al. 2000a), 
NGC 1501 (Sabbadin et al. 2000b, Ragazzoni et al. 2001), 
NGC 6565 (Turatto et al. 2002), NGC 6818 (Benetti et 
al. 2003) and NGC 7009 (Sabbadin et al. 2004). Here we 
present the results for NGC 6741. 

NGC 6741 (PN G033.8-02.6, Acker et al. 1992) is a 
compact (main body~7"x5", halo diametcr~15" ; Curtis 
1918, Schwarz et al. 1992), high-surface brightness, high- 
excitation (class 8, Hyung & Aller 1997) PN with a large 
degree of stratification of the radiation. The powering 
star is very hot (log T*>5.22, Pottasch 1981, Pottasch 
& Preite-Martinez 1983, Heap et al. 1989, Tylenda ct al. 
1989, Kaler & Jacoby 1989) and faint (m v ^19.5, Pottasch 
1981, Tylenda et al. 1989, Kaler & Jacoby 1989, Hyung & 
Aller 1997). 

The [O III] and [N II] apparent morphology of NGC 
6741 (sometimes called "Phantom Streak Nebula") is 
shown in Figs. 1 and 2 (HST-WFPC2 frames retrieved 
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from NASA public archives). Note the vertical (equato- 
rial?) inhomogeneous strip of absorbing knots in the [O 
III] image of Fig. 1, and in Fig. 2 the roundish halo and 
the series of weak, radial [O III] rays punching the [N II] 
skin along and close to the apparent major axis. A multi- 
color HST reproduction of the nebula is given by Hajian 
& Terzian at http://ad.usno.navy.mil/pne/gallery.html 

We immediately point out that the optical appearance 
of NGC 6741, the co-existence of ionic species with a large 
range of ionization potential, IP (from [O I] IP=0 eV to 
[Ne V] IP=97.1 eV; Aller et al. 1985, Hyung & Aller 1997), 
and the characteristics of the central star (hot and faint) 
are suggestive of a recombining PN, i. e. the star has ex- 
hausted the hydrogen-shell nuclear burning and is rapidly 
fading in luminosity (and temperature); the UV flux being 
unable to fully ionize the gas, the outer nebular regions re- 
combine, producing the faint halo (Tylenda 1986, Phillips 
2000). Under many respects NGC 6741 is very much like 
NGC 6565 (Turatto et al. 2002). 

To deepen the kinematics, physical conditions, ionic 
and spatial structure, distance and evolutionary status of 
NGC 6741, we have secured long-slit ESO NTT+EMMI 
echellograms (spectral range AA3900-7900 A, spectral res- 
olution A/AA=R=60 000) at nine position angles (PA). 
The spectra were reduced and analysed using our recon- 
struction 3-D technique (Turatto et al. 2002, Benetti et 
al. 2003). 

The results are presented in this paper, whose plan is 
as follows: Sect. 2 illustrates the observational procedure 
and the reduction method, Sect. 3 defines the kinematics 
of the ionized gas, Sect. 4 quantifies the interstellar and 
circum-nebular absorption, Sect. 5 provides the nebular 
distance, size and age, in Sect. 6 we discuss the parame- 
ters of the central star (temperature, luminosity, mass and 
evolutionary phase) and in Sect. 7 the nebular parameters 
(T c , N c , ionic mass and structure, chemical abundances, 
photo- ionization model), Sect. 8 re-builds the 3-D spatio- 
kinematical structure, Sect. 9 contains the general discus- 
sion, and Sect. 10 draws some conclusions. 

2. Observations and reductions 

NGC 6741 was observed with ESO NTT + EMMI (cchelle 
mode; grating #14, grism #3) at nine equally spaced PA, 
under photometric sky conditions and seeing ranging be- 
tween 0.50" and 0.70". The spectrograph slit (1.0" wide 
and 30" long) was centered on the nebular image, the ex- 
citing star being invisible in the slit-viewer. The echel- 
lograms (exposure time 600s) cover the spectral range 
AA3900-7900 A with resolution R~60 000, and provide the 
kinematical structure of the main ionic species within the 
nebular slices covered by the slit. Bias, zero-order flat field 
and distortion corrections, and wavelength and flux cali- 
brations were performed according to the straightforward 
procedure fully described by Turatto et al. (2002). 

We stress the fundamental difference between our 
frames, covering 80 echelle orders, and the observing pro- 
cedure usually adopted for extended objects, inserting an 



interference filter to isolate a single order. The same is 
valid for the reduction and investigation methods: so far 
long-slit echellograms have been used to obtain, either the 
kinematics in a few ions (in general, [O III], H I and [N 
II]), or the "average" nebular properties (physical condi- 
tions, ionic and chemical abundances) integrated over the 
whole slit length. On the contrary, we proceed to a de- 
tailed, pixel-to-pixel, flux and velocity determination for 
a number of nebular lines, thus inferring the spatial dis- 
tribution of the kinematics, diagnostics, ionic and total 
abundances at the same time. 

The richness of physical information contained in the 
echellograms is illustrated in Figs. 3 and 4, presenting the 
spectral structure of twelve ionic species (from [O I], IP=0 
eV to [Ar V], IP=59.8 eV) at PA=15° and 75° (close to 
the apparent minor and major axes of NGC 6741, respec- 
tively) . 

These figures evidence the spectral characteristics 
common at all PA, in particular: (a) large stratification 
of both the radiation and kinematics (compare, e. g., the 
ionic sequences [O I]-[0 II]- [O III] and [Ar III]-[Ar IV]- [Ar 
V]), and (b) blurred appearance of recombination lines (of 
H I, He I and He II), due to a combination of thermal mo- 
tions, fine-structure and expansion velocity gradient. At 
the same time, they highlight the kinematical and physi- 
cal differences at the two PA: 

- un-tilted, quite uniform emissions close to the appar- 

ent minor axis (Fig. 3), though the blue-shifted gas is 
systematically fainter than the red-shifted gas, 

- tilted, inhomogeneous lines close to the apparent major 

axis (Fig. 4), suggestive of a dense equatorial torus and 
two extended and faint polar lobes. 

In the original NTT+EMMI frames, A7320.121 A of [O 
II] (top- right panel in Figs. 3 and 4; hereafter A7320 A) is 
partially in blend with A7319.044 A, the bluest component 
of the [O II] red quartet (hereafter A7319 A). Since the 
intensity ratio I(A7320 A)/I(A7319 A)=constant=3.071 in 
the N c range here considered (De Robertis et al. 1985, 
Keenan et al. 1999, Sharpee et al. 2004), the A7319 A 
suppression was obtained in the simple fashion illustrated 
in Fig. 5. 

We note in passing that the "relative" spectral and 
spatial resolutions of our echellograms are similar to each 
other: SS=I/ cxp /AV~3 to 5 ~ RR=r/Ar. This means that 
the spatial information of NGC 6741 is as accurate as the 
kinematical information, and both are close to the lower 
limit for the tomographic analysis (Ragazzoni et al. 2001). 
Ergo, the spectral images in Figs. 3 and 4 are affected by 
non-negligible blurring agents. 

Seeing and guiding are the broadening components 
along the spatial (i. e. vertical) axis; in our case the full- 
width at half maximum, W(spat) see ing+guidmg, results to 
be 0.60" for PA=15° (Fig. 3), and 0.65" for PA=75° (Fig. 
4). The true, "intrinsic" profile of the nebula has: 

W(spat)f ntrinslc = W(spat)^ bs - W(spat)^ eeing+guiding . (1) 
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Fig. 3. Detailed spectral image of twelve ionic species in NGC 6741 at PA=15° (close to the apparent minor axis), 
arranged in order of increasing IP (top-left to bottom- right). The original fluxes arc multiplied by the factor given 
in parenthesis, to make each emission comparable with A5007 A of [O III]. The blue-shifted gas is to the left. The 
top-left frame also shows the [O I] night-sky emission at A6300.304 A. In the top-right frame the [O II] line at A7319 
A, partially in blend with A7320 A, has been suppressed (for details, see the text and Fig. 5). 
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Fig. 4. Same as Fig. 3, but for PA=75° (close to the apparent major axis of NGC 6741). 
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Fig. 5. De-blending procedure for the [O II] AA7319-7320 
A doublet (using IRAF packages). Top panel: Image A= 
part of the echcllogram (at PA=75°) centered on A7319 
A. Middle panel: Image B= part of the echellogram (at 
PA=75°) centered on A7320 A. Bottom panel: Image C = 
Image A - (Image B/3.071), where 3.071 corresponds to 
I(7320)/I(7319). 



Concerning the velocity (i. e. horizontal) axis, we must 
take into account: 

(1) instrumental resolution, corresponding to a gaussian 
profile with W(vc1)emmi=5.09 km s^ 1 (measured in 
the [O I] night-sky line at A6300 A); 

(2) thermal motions, generating a gaussian distribution 
with W(vel)thcrmai ^21 .6 x 1(T 2 x T 5 x m" 05 km s" 1 , 
where m is the atomic weight of the element (Clegg et 
al. 1999 and references therein); we use the radial T e 
profile given in Sect. 7; 

(3) turbulence, i.e. random, small-scale motions; this is 
a very uncertain parameter for PNe, whose quantifi- 
cation is deferred to a dedicated paper (in prepara- 
tion), based on very-high resolution (R~115 000) spec- 
tra secured with Telescopio Nazionale Galileo (TNG) 
+ SARG on a representative sample of targets. In the 



Fig. 6. High-contrast [O III] and [N II] line profiles (in log- 
arithmic scale) of NGC 6741 close to the apparent minor 
and major axes (PA=15° and 75°, respectively) showing 
the faint halo and red-shifted tail. Same orientation as 
Figs. 3 and 4. 



present case of NGC 6741, at the moment we can only 
infer that W(vel) tU rb is below 10.0 km s _1 , as indi- 
cated by the sharpness of the spectral images of for- 
bidden lines along the velocity (i. e. x) axis (see Figs. 3 
and 4; larger turbulences should produce blurred spec- 
tral images of forbidden lines, similar to the Ha one). 
According to the "a posteriori" analysis presented at 
the end of Sect. 3, in the following we will assume 
W(vel) t urb= 3.5 km s _1 , in partial agreement with the 
general results by Neiner et al. (2000), and Gesicki et 
al. (2003) (but see the caveat in Sect. 3); 
(4) fine-structure (only for recombination lines); follow- 
ing Clegg et al. (1999), W(vel) fino _ s . -7.5 km s" 1 for 
Ha; moreover, we adopt W(vel)fi nc - S .=5.0 km s _1 for 
A4686 A (after suppression of the 20 km s 1 blue- 
shifted tail with the method illustrated in Fig. 5) and 
A6560 A of He II, and A5876 A of He I. 

Note that thermal motions and/or fine-structure repre- 
sent the main broadening factors for recombination lines, 
whereas instrumental resolution, thermal motions and tur- 
bulence are (more or less) equivalent for forbidden lines. 
The final broadening along the velocity axis is: 



W(vcl)L r =E W ( vcl )i 



i=l 



and the true, intrinsic full- width at half maximum: 



W(vel) 



intrinsic 



W(vel)L-W(vel) 



blur' 



(2) 



(3) 
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W(vel)i ntr insic (in krn s -1 ) corresponds to the expansion 
velocity range of the emitting layer; it can be trans- 
formed into arcsec by means of the general expansion 
law of the ionized gas: in NGC 6741 being T4 xp (km 
s _1 )=13(±l)xR" (see Sect. 3), we have W(vel) intr insic 
(arcsec)= [W(vel) intr insic (km s _1 )]/13. 

Let us consider, for example, the spectral images at 
PA=15° (Fig. 3). The intrinsic FWHM in both the zvpc 
and cspl (as defined at the opening of Sect. 3) results to 
be 0.52" to 0.72" (forbidden lines), 0.77" to 0.95" (He 
I), 1.00" to 1.15" (He II), and 1.35" to 1.45" (H I), thus 
confirming that: 

(a) each emitting region is extremely sharp, 

(b) the whole spectral images must be carefully decon- 
volved: to this end we use the Richardson-Lucy algorithm 
(Richardson 1972, Lucy 1974), and a point-spread func- 
tion given by a bi-dimensional gaussian profile character- 
ized by W(spat) se eing+guiding and W(vel) b i U r- 

So far we have considered the bright main shell of NGC 
6741. At lower-flux cuts the faint halo appears (Fig. 6), 
whose spectral signature supports the recombination hy- 
pothesis (note, in particular, the broad, un-tilted halo- 
emission at PA=15°, and the broad, tilted halo-emission 
mimicing the kinematics of the main nebula at PA=75°). 
In fact, according to the current evolutionary models (e. 
g. Schonberner et al. 1997, Steffen ct al. 1998, Marigo et 
al. 2001, Pcrinotto et al. 2004a), a PN halo can be: 

- an AGB halo, i. e. the low-density envelope of an op- 

tically thin PN, directly ionized by the UV flux of the 
bright central star. It represents the AGB mass-loss in 
the prc-supcrwind phase, and expands at the original 
ejection velocity (V cxp ~10-15 km s _1 ), independent 
on the kinematical field of the main nebula; 

- a recombination halo, corresponding to the outer, al- 

most neutral parts of the main nebula, no more reached 
by the UV flux of a hot central star which has ex- 
hausted the H-shell nuclear burning and passed the 
turnover point in stellar evolution, rapidly fading in 
luminosity and temperature. The recombining layers 
essentially retain the general kinematical properties of 
the main thick nebula (as observed in the halo of NGC 
6741). 

Fig. 6 also shows the faint red-shifted tail present in 
the strongest emissions of NGC 6741 at all PA, whose na- 
ture (halo? interaction with the ambient ISM? early- AGB 
wind? instrumental scattered light?) remains unclear. 

This preliminary, qualitative description enhances the 
complexity of NGC 6741 (a common characteristic of all 
PNe analysed at adequate spatial and spectral resolu- 
tions). Before starting a thorough investigation of the 
kinematical and physical information contained in the 
spectra, we underline that the general characteristics of 
NGC 6741 - in particular: (a) nebular compactness, and 
(b) star weakness (the stellar continuum, generally used 
as a position marker, is absent in the cchcllograms) - do 
make this nebula the ideal target to test the reliability 
limits of our 3-D reconstruction procedure. 



3. The gas spatio-kinematics 

According to Sabbadin et al. (2004, and references 
therein), the overall spatio-kinematical properties of a reg- 
ularly expanding nebula can be derived by combining the 
kinematics of the gas projected at the apparent position 
of the star (the "central star pixel line" , cspl, of the echel- 
lograms, common at all PA) with the spatial profile at the 
systemic radial velocity (the "zero velocity pixel column" , 
zvpc, at each PA). 

Though the stellar continuum is absent in the cchcl- 
lograms of NGC 6741, we obtain a satisfactory cspl- 
location and zvpc-alignement (within ±0.1") of the echelle 
orders thanks to: (a) the central, symmetrical position of 
the star in the broad-band HST-WFPC2 nebular images, 
(b) the cspl- and zvpc-calibration given by the continuum 
spectrum of the central star of NGC 7009 (observed at 
the same night, and with the same instrumental configu- 
ration). 

The peak separations in the cspl of NGC 6741, 2V cxp , 
are contained in Table 1 (last column), where ions are 
put in order of increasing IP. Typical errors are 1.0 km 
s _1 for the strongest forbidden emissions (like A4959-5007 
A of [O III] and A6548-6584 A of [N II]) to 2.0 km s" 1 
for the faintest ones (in particular, A6312 A of [S III], 
A471 1-4740 A of [Ar IV] and A6435-7005 A of [Ar V]). 
The uncertainties for recombination lines are: 2.0 km s _1 
for A4861-6563 A of H I, A4686-6560 A of He II and A4640 
A of N III, and 1.5 km s" 1 for A5876-6678 A of He I. 

The agreement with the kinematical results given in 
the literature (also contained in Table 1) is quite poor, 
probably due to the compactness of NGC 6741, combined 
with image rotation on the slit (long-exposure Coude spec- 
tra by Wilson 1950) or inadequate spatial resolution (cir- 
cular aperture 18" in diameter by Robinson et al. 1982). 

Concerning the zvpc at the nine PA of NGC 6741, the 
intensity peak separations, 2r zvpc , in the different ionic 
species are listed in Table 2. 

To assemble the kinematical results of Table 1 and the 
spatial results of Table 2 we can follow two ways: 

(I) Search for the PA with R zvpc ~ R csp i. Let us as- 
sume for the main shell of NGC 6741 the most gen- 
eral spatial structure, i. e. a tri-axial ellipsoid with 
axes a, b and c. At PA=75° to 115° (close to the ap- 
parent major axis, Fig. 4) the line-tilt indicates that 
the observer is not aligned with the major axis, and 
the overall emission structure (suggestive of a dense 
equatorial torus + two faint polar lobes) excludes the 
oblate (a=b) ellipsoid hypothesis. Moreover, the ab- 
sence of line-tilt at PA=15° (perpendicular to the ap- 
parent major axis, Fig. 3) means that either the minor 
axis of the ellipsoid is close to PA=15°, or the inter- 
mediate axis of the ellipsoid is close to PA=15°, or the 
ellipsoid is prolate (b=c). In all cases we conclude that, 
in first approximation, R(minor axis)~R(intermcdiatc 
axis)~R zvpc (PA=15°)~ R CS pi, providing (through 
Tables 1 and 2): V cxp (km S - 1 ) = 13(±2) xR". 
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Table 1. Peak separation in the cspl of NGC 6741 



Ion 


IP range (eV) 




2K xp (km/s) 




Wilson (1950) 


Robinson et al. (1982) 


this paper 


[0 I] 


0.0-13.6 






55.1 


[S II] 


10.4-23.3 






54.0 


[0 II] 


13.6-35.1 


44.2 




53.7 


H I 


>13.6 






44.0 


[N III 

L J 


14.5-29.6 


42.1 




53.5 


[s ml 

IV J 


23.4-34.8 






48.2 


He I 


24.6-54.4 






48.0 


[Ar III] 


27.6-40.7 






47.6 


[O III] 


35.1-54.9 


41.6 


42 


46.1 


[Ar IV] 


40.7-59.8 






33.7 


[Ne III] 


41.0-63.4 


41.0 




48.2 


N III 


47.4-77.5 






40.0 


He II 


>54.4 






34.0 


[Ar V] 


59.8-75.0 






25.0 


[Ne V] 


97.1-126.2 


0.0: 







[O III] [O III] [O III] 
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- lnMi + 



Fig. 7. Nebular tomography vs. spectral image connection. In this example we use the [O III] spectral images of NGC 
6741 at PA=75° and PA=15° (same orientation as Figs. 3 and 4). Left panel: tomographic map of the nebular slice 
intercepted by a spectrograph slit aligned with the apparent major axis of an ellipsoidal PN denser at the equator 
than at the poles; a and b are the polar and equatorial axes, respectively, of the elliptical nebular slice. Central panel: 
spectral image of the nebular slice shown in the left panel; a' and b' are the polar and equatorial axes, respectively, of 
the spectral image. Right panel: [O III] tomographic reconstruction of NGC 6741 at PA=15°. See the text for details. 



(II) Spectral image deformation along the major axis. 

Let us consider a generic, expanding (I4 X p ocR) tri- 
axial ellipsoid denser in the equatorial region than 
at the poles, seen at an intermediate direction. A 
spectrograph slit aligned with the apparent major 
axis intercepts the radial slice of nebula shown in 
Fig. 7 (left panel): an ellipse with polar axis a and 



equatorial axis b, perpendicular to each other. The 
corresponding spectral image is an ellipse too (Fig. 
7, central panel), but deformed (with respect to the 
tomographic map) by the telescope + spectrograph 
characteristics. Note that: 
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Ion 2r zvpc (arcsec) 



PA=15° 35° 55° 75° 95° 115° 135° 155° 175° 



fO 11 

J.J 


4.3 


4.6 


5.7 


7.3 


[S II] 


4.1 


4.5 


5.5 


7.1 


[0 II] 


4.0 


4.3 


5.2 


6.6 


HI 


3.5 


3.5 


4.1 


6.0 


[N II] 


4.1 


4.4 


5.4 


6.8 


[S III] 


3.6 


4.1 


4.6 


6.0 


He I 


3.8 


4.0 


4.5 


6.2 


[Ar III] 


3.7 


3.8 


4.5 


5.9 


[0 III] 


3.6 


3.7 


4.4 


5.3 


[Ar IV] 


2.7 


3.2 


3.6 


3.9 


[Ne III] 


3.6 


3.5: 


4.0 


4.5 


N III 


3.0: 


3.1: 


3.6: 


4.0: 


He II 


2.6 


2.9 


3.4 


3.8 


[Ar V] 


2.3 


2.6 


2.8 


3.0 



6.8 


7.1 


6.1 


5.6 


4.6 


fi fi 




u.u 


U.O 


4.4 


6.5 


6.6 


5.8 


5.2 


4.3 


6.0 


6.0 


5.2 


4.5 


3.8 


6.5 


6.5 


5.9 


5.2 


4.4 


5.7 


6.0 


5.5 


4.6 


3.9 


5.7 


6.0 


5.5 


4.6 


4.0 


5.0 


6.0 


5.5 


4.6 


4.1 


4.7 


5.5 


5.0 


4.4 


3.8 


3.9 


5.3 


4.5: 


3.4 


3.1 


5.0 


5.6 


5.0 


4.0 


3.6 


4.2 


5.1 


4.4: 


3.4: 


3.2: 


3.9 


5.0 


3.9 


3.2 


2.9 


3.5 


4.0 


3.6 


3.0 


2.7 



- the polar axis (a') and the equatorial axis (b') of 
the spectral image are no more perpendicular to each 
other, 

- the original tomographic map (Fig. 7, left panel) can 
be obtained by means of a simple compression of the 
spectral image (Fig. 7, central panel) along the x (i. e. 
velocity) axis, till having a'_Lb'. 1 

In practice, and reversing the foregoing procedure, we 
compress the [O III] spectral image of NGC 6741 at 
PA=75° (Fig. 7, central panel) by a factor 1.60 along 
the x axis, thus correcting for the spectral deforma- 
tion introduced by NTT+EMMI on the nebular slice 
intercepted by the slit, and obtaining: 

- the tomographic map shown in Fig. 7 (left panel), 

- the expansion law T4 X p(km s _1 )=13(±l) xR". 
This is the same law given by method (I). In fact, the 
compression factor along the x axis being the same at 
all PA, we can repeat the procedure for the [O III] 
spectral image of NGC 6741 at PA=15° (close to the 
apparent minor axis; Fig. 3), thus recovering the to- 
mographic map shown in the right panel of Fig. 7, i. 
e. an almost circular ring with R zvpc (PA=15°)~ R CS pi, 
QED! 2 

The detailed cspl-zvpc relation for NGC 6741 at PA=15°, 
shown in Fig. 8, indicates that: 

a) on the whole, the ionized gas follows Wilson's law: 
the high-excitation zones expand more slowly than the 
low-excitation ones, and a positive correlation exists be- 
tween the expansion velocity and the size of the monochro- 
matic image (Wilson 1950); 

1 This visual effect can be seen by rotating Fig. 7 around the 
y axis. 

2 quod erat demostrandum; Latin for: which was to be 
proved. 



b) the range of both 2r zvpc and 2V CKp is quite large 
(2.3 to 4.3 arcsec and 25.0 to 55.1 km s _1 , respectively), 
suggesting a broad radial density profile and large strati- 
fication of the radiation; 

c) the general expansion law, V cxp (km s _1 )=13xR", 
fails in the innermost, highest ionization layers marked by 
the [Ar V] emissions: they expand slower than expected. 

Point c) is quite peculiar among PNe, and deserves 
some comments. Deceleration is present (and clear) in 
both the A6435 A and A7005 A [Ar V] lines at all PA of 
NGC 6741 (except along and close to the apparent major 
axis, where the open-ended structure of the [Ar V] spec- 
tral image prevents any conclusion; see Fig. 4, lower-right 
panel), whereas no evidence of deceleration appears in 
other high ionization species, like He II (due to the blurred 
emissions) and N III (A4640 A is too faint). Though a de- 
tailed spatio-kinematical study at even higher ionization 
stages appears indispensable - for example in the forbid- 
den line of Nc+ 4 (IP=97.1 cV) at A3425 A (which is out- 
side our spectral range) - a support to the deceleration hy- 
pothesis comes from the classical work by Wilson (1950), 
who obtained 14x P [Ne V]~0 km s _1 (quite uncertain). 

We believe NGC 6741 is the first PN showing clear evi- 
dences of deceleration in the highest ionization layers (the 
second candidate being IC 418: a preliminary analysis of 
our NTT+EMMI echellograms, taken at six PA, suggests 
the probable presence of inf ailing gas). 

Vice versa, recent results by Gesicki & Zijlstra (2003, 
3 PNe) and Gesicki ct al. (2003, 14 PNe) - based on high- 
resolution emission profiles integrated along the slit, and 
spherical shell hypothesis for the emitting gas - indicate 
that acceleration is a quite common property of the PN 
innermost layers (i. e. "U" -shaped expansion profile), due 
to the dynamical contribution by the shocked, hot wind 
from the central star. Gesicki & Zijlstra (2003) adopted 
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A4686 A of He II as diagnostic of the high-ionization kine- 
matics, whereas Gesicki et al. (2003) used A5007 A of [O 

III] (in a few cases A6560 A of He II). 
The following caveat are in order: 

1) the spherical shell assumption is wide of the mark; 

2) A5007 A of [O III] (IP range 35.1 to 54.9 eV) is a 
poor diagnostic of the highest ionization strata (except for 
very-low excitation PNe); 

3) the recombination lines of hydrogen and helium 
do suffer severe blurring effects (thermal motions, fine- 
structure and expansion velocity gradient across the neb- 
ula combined with the small number of ionization stages 
of H and He), introducing spurious kinematical results. 

In particular, A4686 A, He II Paschen a, con- 
sists of thirteen fine-structure components spread in the 
AA4685. 377-4685.918 A range (corresponding to AV=34.6 
km s _1 ), with a strong blue-shifted tail (see Figs. 3 
and 4). In the case of A6560 A, He II Pickering j3, 
the fine-structure components are even nineteen, cover- 
ing the AA6559.769-6560.209 A spectral range (AV=20.1 
km s _1 ). 

To further test this point, let us consider the emis- 
sion line profiles integrated along the slit of NGC 6741 at 
PA=15°. The results for A4686 A and A6560 A of He II, 
and the ionic sequence of argon - [Ar III] at A7135 A, [Ar 

IV] at A4740 A and [Ar V] at A7005 A - are shown in Fig. 
9; they confirm that: 

a) spherical symmetry is a simplistic assumption, 

b) the He II recombination lines fail to fit the well- 
defined stratification of both the radiation and expan- 
sion present in the ionic sequence of argon, FWHM(Hc 
II, IP>54.4 eV) being intermediate between FWHM([Ar 
III], IP range 27.6-40.7 eV) and FWHM([Ar IV], IP range 
40.7-59.8 eV), whereas, according to the detailed kine- 
matical results presented in this section, we would expect 
FWHM(He II) < FWHMQAr IV]) « FWHMQAr III]). 

The same discrepancy occurs when considering either 
more parameters (like full- width at 10% maximum flux), 
or different ionic sequences, or further PA of NGC 6741, 
or other PNe of our sample. 

A direct confirmation of the misleading kinematical 
results - in particular, the high-velocity of the innermost 
layers - provided by the combination (spectral profile of 
recombination lines integrated along the slit) + (spherical 
symmetry assumption) comes from a comparative analysis 
of the emission line profiles contained in Gesicki & Zijlstra 
(2003) and Gesicki et al. (2003): in all cases (17 PNe) the 
forbidden lines of the ionic species with the highest IP, 
also show the sharpest emission profile (as expected for a 
simple, positive V CKp vs. radius relation). 

All this questions the general validity of the spatio- 
kincmatical studies based on the spherical symmetry as- 
sumption and/or spectral profiles of recombination lines 
integrated along the slit, and weakens (or even cancels) 
the reliability of their results on turbulence, radial kine- 
matics, matter distribution and ionization, whose quantifi- 
cation needs detailed studies at high spatial and spectral 
resolutions. 



m 
o 



> 
cv? 




2 4 
2r(zvpc) [arcsec] 



Fig. 8. The cspl-zvpc relation for NGC 6741 at PA=15°, 
superimposed to the adopted expansion law V exp (km s _1 ) 
=13 xR", which is valid across the whole nebula, except 
in the innermost, highest-ionization, decelerated regions 
marked by [Ar V]. 



Let us go on. Fig. 10 (multi-color map only in the 
"free" electronic version of the paper, for lack of funds) 
shows the complete velocity field at the nine observed 
PA of NGC 6741. We select He II, [O III] and [N II] 
as markers of the high, medium and low-excitation re- 
gions, respectively. These position-velocity (P-V) maps 
are relative to the systemic heliocentric velocity of the neb- 
ula, Vrad0= +39.6(±1.0) km s _1 , corresponding to Vlsr= 
+56.4(±1.0) km s _1 , and are scaled (along the x axis) ac- 
cording to V exp (km s _1 )~13xR", i.e. they reproduce the 
tomographic maps of the nebular slices covered by the 
spectrograph slit. 

Fig. 10 highlights the kinematical complexity and 
large stratification of the radiation within NGC 6741. 
The main nebula consists of an almost-prolate ellipsoid 
(a~7.4", a/b~1.8, a/c~2.0), whose major axis (projected 
at PA~95°) forms an angle of 55(±3)° with the line of 
sight: on the whole, the eastern part of the nebula is ap- 
proaching the observer, and the western part receding. 

All this throws new light on two specific fields: (I) na- 
ture of the curious skip of absorbing knots present in Fig. 
1 (upper panel), and (II) turbulent motions. 

(I) skip of absorbing knots: the central location, 
orientation (vertical), and light curvature (concavity to- 
wards East) are suggestive of an inhomogeneous belt of 
neutral matter (gas + molecules + dust) embedding the 
dense ionized gas of the equatorial regions. The amount 
of circum-nebular neutral gas can be estimated from the 
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Fig. 9. Selected emission line profiles integrated along the 
slit for NGC 6741 at PA=15°, showing the spurious kine- 
matical results provided by the recombination lines of He 
II. Symbols: continuous linc= [Ar III] at A7135 A; dot- 
ted line= [Ar IV] at A4740 A; short-dashed line= [Ar V] 
at A7005 A; long-dashed linc= He II at A4686 A; dotted- 
dashed line= He II at A6560 A. For details, see the text. 



[O III] flux-depletion suffered by the underlying ionized 
nebula, being: 

F(5007) off _knot , . , mm 

l0 S pfmn?^ = fc 5007 - KH 3 +c(H^) circum _ ne bular, (4) 

V (OUUf Jon-knot 

where k\ is the extinction coefficient (Seaton 1979), 
and c(H/?) c i rcum _ no buiar the logarithmic extinction at H/3 
caused by the local absorbing matter. 

Intensity scans in the neighbourhood of the deep- 
est knots provide p^oo?)"":^ U P to 1A (±0-1), i- e. 
c(H/3) circum _ nebular up to 0.18 (±0.03). Assuming a "nor- 



mal" gas-to-dust ratio (Spitzcr 1978, Bohlin et al. 1978) 
and c(H/?)=1.48xE(B-V) (Acker 1978), we have: 



n(HI) 



4.8 x 10 21 x E(B - V) 



circum— nebular 



Al 



(5) 



where n(HI) is the H I density (atoms cm -3 ), and 
Al the radial thickness (cm) of the absorbing layer. For 
D=2000 pc (see Sect. 5) and Al=[r(halo) - r(main nebula 
at PA=15°)]/2~2.0" (from Figs. 1 and 2), we obtain an in- 
dicative value of n(HI) ~7xl0 3 atoms cm~ 3 . Such a high 
density of the circum-nebular matter in NGC 6741 is a 
further support to the recombination hypothesis outlined 
in the previous sections. 



(II) turbulent motions. In Fig. 10 the tomographic 
map close to the apparent minor axis (i.e. at PA=15°) be- 
ing almost circular and quite homogeneous, we can assume 
W(spat)intrinsic(zvpc at PA=15°) ~ W(vel) intr insic(cspl), 
thus inferring W(vel) tU rb (through Eqs. (1) to (3), Vo Xp (km 
s _1 ) = 13xR", and the T c radial profile given in Sect. 
7.1). We overlook the recombination lines of hydrogen 
and helium (dominated by thermal motions and/or fine- 
structure), and consider the strongest forbidden lines. The 
analysis of A6300 A ([O I]), A6731 A ([S II]), A6584 A ([N 
II]), A7135 A ([Ar III]), and A5007 A ([O III]) provides 
W(vcl)turb= 3.5(±2.0) km s _1 , with no evident relation to 
the ionization degree. Thus, in spite of the crude assump- 
tions and wide uncertainties, we conclude that turbulent 
motions in NGC 6741 are quite modest. 

4. The absorption (interstellar + circum-nebular) 

In general, the observed line intensities must be corrected 
for absorption according to: 



log^2IL = fc A x c (H/3) tot , 



1(A) 



obs 



(6) 



where k\ is the extinction coefficient (Seaton 1979), and 
c(H/3) to t the logarithmic extinction at H/3 given by: 



c(H/3) to t — c(H/3)i ntorst cllar + c(H/3) circum 

— nebular 



(7) 



with c(H/?) circum _ nebu i ar =0 for an optically thin, fully 
ionized, density bounded nebula. Decidedly, this is not 
the case of NGC 6741, which is optically thick, ionization 
bounded, and wrapped up in a dense cocoon of almost- 
neutral matter. 

The extinction estimates reported in the literature 
(mean values along the nebular slice covered by the spec- 
trograph slit) cluster around c(H/3) to t=l-05 (Kaler & Lutz 
1985, Kaler & Jacoby 1989, Cahn et al. 1992, Hyung & 
Aller 1997). 

To disentangle the complex absorption over NGC 6741, 
we apply the Ha/H/3 analysis of the whole spectral image, 
as introduced by Turatto et al. (2002). Fig. 11 shows the 
F(Ha)/F(H/3) isophotal contours superimposed to the H/3 
spectral image for three representative PA of NGC 6741: 
at PA=15° (close to the apparent minor axis; untiltcd H/3 
spectral image) F(Ha)/F(H/3) increases outwards along 
the spatial (i. e. y) axis, peaking beyond the top of the 
H/3 flux. At PA=55° (intermediate PA; tilted H/3 spectral 
image) the F(Ha)/F(H/3) raise along the y axis is over- 
whelmed by the broad maximum at S-W of the nebular 
centre (at the expected position of the absorbing belt vis- 
ible in Fig. 1, upper panel). The same occurs at PA=75° 
(close to the apparent major axis; tilted H/3 spectral im- 
age), where F(Ha)/F(H/3) peaks at W-SW of the centre, 
in correspondence of the equatorial absorbing belt. 

Summing up the Ha/H/3 intensity distribution at the 
nine observed PA of NGC 6741, we infer that: 

a) c(H/3) 

interstellar 

0.95 (±0.05), 
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Fig. 10. Combined position-velocity maps in the nine observed PA of NGC 6741 at high (He II), medium ([O III]) 
and low ([N II]) excitation (multi-color maps in the electronic version of the paper; blue=He II, green=[0 III], and 
red=[N II]), scaled according to the relation V cxp (km s _1 ) = 13xR". The orientation of these tomographic maps is the 
same of Figs. 3 and 4. 



b) c(H/3) circum _ ncbular is 0.10 (±0.03) in the central region 
of the nebular image (out of the equatorial absorbing 
belt), and raises up to 0.20 (±0.05) within the equa- 
torial absorbing belt, and at the edge of the ionized 
zone. 



Moreover, a modest decrease of F(Ha)/F(H/3) in the in- 
nermost regions (out of the equatorial absorbing belt) is 
suggestive of a local increase of T c . 



Last, an indicative value of the total neutral mass 
embedding the ionized nebula (through Eq. (5) and us- 
ing simple geometrical considerations) is M nou t ra i— 0.20 
(±0.05) M Q . 

Though the small angular size prevents a deeper anal- 
ysis, these results do represent the umpteenth sign of 
the complex structure and peculiar evolutionary phase of 
NGC 6741, and will be a precious support in the determi- 
nation of the topical parameter, i. e. distance. 
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Fig. 11. F(Ha)/F(H/3) isophotal contours superimposed to the H/3 spectral image for three representative PA of NGC 
6741. Left panel: PA=15° (close to the apparent minor axis); central panel: PA=55° (intermediate direction); right 
panel: PA=75° (close to the apparent major axis). Same orientation as Figs. 3 and 4. The isophotal contours cover 
the range 5.20 (the outermost) to 7.60, with a constant step of 0.30. 



5. The nebular distance, size and age 

The "statistical" distance of NGC 6741, provided by two 
dozen catalogues using different methods and assump- 
tions, is: 

<D>(Shklovsky)~3300(±1000) pc 

<D>(ionized mass-radius relation)~1500(±700) pc 

<D> (other mcthods)~2000(±1000) pc, 

where <D>(Shklovsky) represents an upper limit of the 
true distance, NGC 6741 being an optically thick, ioniza- 
tion bounded PN. 

Individual values reported in the literature are: 1300 pc 
(nebular radial velocity combined with the circular law of 
galactic rotation; Acker 1978), 2100 pc (Acker 1978) and 
1400 pc (Pottasch 1983) (both based on the large-scale 
galactic interstellar absorption map by Lucke 1978), and 
1500 pc (color-excess vs distance relation for early- type 
stars within 1.5° of the nebula; Kaler & Lutz 1985). 

We tried to determine the nebular parallax by combin- 
ing the expansion velocity field with the angular expan- 
sion measured in first- and second-epoch HST-WFPC2 
frames. The NASA public archives contain 21 images of 
NGC 6741, taken at two epochs separated by 2.97 years 
(programs GO 7501 and GO 8773; P. I. Arsen Hajian). 
The [O III] and [N II] multiple exposures (WFPC2 cen- 
tral planetary camera, PC; pixel sizc=0.0455 arcsec) were 
co-added, corrected for optical distortions, aligned and ro- 
tated using IRAF packages (Reed ct al. 1999; Palen ct al. 
2002). 



No apparent image-shift is detected, i. e. imagc- 
shift<ipixel, and ^|<8xl0" 3 arcsec yr _1 . For an opti- 
cally thin, fully ionized PN this implies D>1400 pc, being: 



D(pc) = 



0.211T/ cxp (km s' 1 ) 
H (arcsec yr _1 ) 



(8) 



But NGC 6741 is optically thick (probably recombining), 
and Eq. (8) cannot be applied, since in this case we should 
compare the expansion velocity of the ionized gas with 
the angular expansion of the ionization front. To be no- 
ticed that, in the limit case of a nearby nebula in a very 
deep recombination phase, we even expect a detectable 
contraction of the ionization front. Thus, we infer (a mere 
speculation, at the moment) that D(NGC 6741)>1500 pc, 
and/or the nebula is at the end of the recombination phase 
(or at the reionization start). 

As a last chance for the nebular distance determina- 
tion, we decided to enrich (and sharpen, if possible) the 
broad color-excess vs. D relation given by Kaler & Lutz 
(1985) with the recent bibliographic reports. The scan- 
ning of the NGC 6741 neighbourhood with the SIMBAD 
facilities of CDS (Centre de Donnees astronomiques de 
Strasbourg) gave fruitful results; besides two dozen (low- 
weight) field stars with accurate photometry and spectral 
type (in most cases the luminosity class is absent; we as- 
sume a luminosity class V), we have identified two impor- 
tant distance markers: 

- the open cluster OCL 88 (at an apparent distance 
of 15.4' from NGC 6741), characterized by E(B-V)=1.0 
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K and A^ c =2xl0 3 cm -3 ), i. e. 3.2<A V <4.7, and, using the 
oscuration law by Forbes (1985), D> 6000 pc for all four 
nebulae. 

The last result appearing quite improbable, and, in 
our opinion, questionable, in the following we will adopt 
the color-excess vs. distance relation shown in Fig. 12. It 
provides D(NGC 6741) = 2000 (±300) pc. 

The linear size of the main nebula is 0.036 pc x 0.020 
pc x 0.018 pc (major, intermediate and minor semi-axes, 
respectively), whereas both the spherical recombining halo 
and the [O III] rays punching the [N II] skin close to the 
apparent major axis extend up to 0.080 pc. 

The "kinematical" age of NGC 6741, t kin =R/V; x p, 
is about 750 years, and the "true" age ttruc— 2R/ 
[V cxp (today)+K X p(AGB)]~ 1400 years. 

Summing up: NGC 6741 is a very young, compact, 
high surface brightness (i.e. dense) PN. The combination 
of: (a) high excitation (up to [Ne V]) of the innermost lay- 
ers, (b) low ionization skin, and (c) almost neutral, large 
density halo, is indicative of a very hot central star at low 
luminosity. A deep analysis of the powering engine is in 
order. 



Fig. 12. Color-excess vs. distance relation in the direc- 
tion of NGC 6741. Symbols: small dots= field stars; 
filled square^ open cluster OCL 88; large dot= classical 
Cepheid V336 Aql. 

and D=3000 pc (Sagar & Griffiths 1998, and referenced 
therein), 

- the classical Cepheid V336 Aql (apparent dis- 
tance from NGC 6741=37.0'), with log P=0.8636 days, 
<m v >=9.875, E(B-V)=0.61-0.64, and D=1800-2100 pc 
(Pel 1976, Feast & Whitelock 1997, Metzger et al. 1998, 
Bcrdnikov et al. 2000). 

The improved E(B-V) vs. D relation in the direction of 
NGC 6741 is shown in Fig. 12: the interstellar absorption 
starts at D~ 200-300 pc, and gradually increases up to 
3000 pc (and beyond). This is in partial disagreement with 
the literature reports. We recall that NGC 6741 is within 
the Aquila-Scutum cloud, a galactic region characterized 
by a high, quite uniform star density, suggestive of a low 
oscuration. According to Forbes (1985), it shows an almost 
immediate climb at Ay=2.3 mag. within 300 pc; beyond 
this distance, there is little or no increase of extinction out 
to D~6000 pc. 

To test this point, let us consider the four more PNe 
projected within the Aquila-Scutum cloud (at an apparent 
distance from NGC 6741 lower than 90'). They are CBSS 
3, K 3-19, M 1-66 and K 3-20 (we reject a fifth object, Sp 
2-151, because of the unreliable line intensities - in par- 
ticular F(Ha)/F(H/3)~1.0 - and peculiar spectral-type A 
for the exciting star reported by Acker et al. 1992, prob- 
ably due to the partial blending with a relatively bright, 
mR~13.0, field star). The four selected PNe in the Aquila- 
Scutum cloud do have 1.5<c(H/3)<2.2 (line intensities by 
Acker ct al. 1992, and Cappcllaro et al. 1994, for T c =10 4 



6. The central star parameters 

The HST-WFPC2 frames of NGC 6741 taken through 
the broad-band filter F555W (central wavelength=5407 
A, bandwidth=1236 A) provide m v =20.09 (±0.05), where 
the un-known star color is the main source of inaccuracy. 
Previous ground-based estimates reported in the litera- 
ture are: 19.5 (Pottasch 1981), 17.6 (Tylenda et al. 1989), 
19.16 (Kaler & Jacoby 1989) and 19.26 (Hyung & Aller 
1997) for m v , and >20.2 (Gathicr & Pottasch 1988) and 
18.2 (Tylenda et al. 1989) for m B . 

The H I and He II Zanstra temperatures are given 
by the stellar magnitude, combined with both the total 
H/3 nebular flux, log F(H/3) obs =-11.32 (±0.03) mWxm- 2 
(Kaler & Lutz 1985, Kaler & Jacoby 1989, Acker et al. 
1991, Cahn et al. 1992, this paper), and the flux ratio 
F(A4686 A)/F(H/?)=0.40(±0.03) (Aller at al. 1985, Kaler 
& Jacoby 1989, Hyung & Aller 1997, this paper). 

We obtain log(T z H I)=5.33(±0.07) and log(T z He 
II)=5.23(±0.07), thus confirming the peculiarity already 
reported by Heap et al. (1989) and Hyung & Aller (1997), 
i. e. the Zanstra discrepancy is reversed. T Z H I>T z He 
II is a typical signature of recombining PNe, e. g. NGC 
6565 (Turatto et al. 2002), where the ionization and ther- 
mal structure are out of equilibrium, and the recombina- 
tion processes dominate. These are faster for the higher- 
ionization species (Tylenda 1986, Stasinska 1989, Marten 
& Szczerba 1997); thus, in the following, we will adopt 
T„~T z He IKT Z H I. 

The stellar luminosity (using D=2000 pc, and the bolo- 
metric corrections by Schonbcrncr 1981) results to be log 
L*/L =2.75(±O.15). 

The high temperature and low luminosity of the star, 
added to the short nebular age, suggest that the stellar 
mass, M*, is larger than the average value (~0.60 M Q ) of 
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the PNe nuclei. According to the evolutionary tracks by 
Schonberner (1981, 1983), Iben (1984), Wood & Faulkner 
(1986), Blocker & Schonberner (1990), Vassiliadis & Wood 
(1994) and Blocker (1995), the 0.66-0.68 M Q post-AGB 
star of NGC 6741 has recently (a few hundreds years ago) 
exhausted the hydrogen-shell nuclear burning, and is fad- 
ing along the white dwarf cooling sequence. 

To be noticed: the early luminosity decline of the 
star was very fast (and caused the nebular recombina- 
tion), but later it gradually slowed, so that wc cannot 
exclude that, at present, the gas has reached (or even 
passed) the contraction-expansion equilibrium condition, 
(l/3)x(d(lnL*)/dt)=-2/tkin in the classical Stromgren 
model, thanks to the matter dilution due to expansion. 

New input will come from the nebular physical con- 
ditions, radial ionization structure and photo- ionization 
model. 



7. The nebular parameters 

7.1. Physical conditions and ionized mass 

Following Sabbadin et al. (2004, and references therein), 
the radial profile of the physical conditions is given by 
the zvpc, which is independent on the expansion velocity 
field, since it represents the gas in the plane of the sky, 
whose motion is tangential. For T c we use the classical 
diagnostic line ratios of ions in p 2 and p 4 configurations 
(A5007 A/A4363 A of [O III] and A6584 A/A5755 A of 
[N II]; Aller 1984, Osterbrock 1989), whereas iV comes 
from both diagnostic line ratios of ions in p 3 configuration 
(A6717 A/A6731 A of [S II] and A4711 A/A4740 A of [Ar 
IV]), and the absolute Ha flux distribution. 

The main limitation is connected to the NGC 6741 
compactness, large stratification of the radiation, and 
weakness of the [N II] auroral line and the [S II] and 
[Ar IV] doublets: T C [N II], N C [S II] and V c [Ar IV] can 
be derived only at (or close to) the intensity peak of the 
corresponding emission. 

Moreover, the large Ha broadening implies a com- 
plex deconvolution for instrumental resolution plus ther- 
mal motions plus fine-structure (see Sect. 2), lowering the 
accuracy of the F(Ha) zvpc and N c (Ha) profiles. Thus, ac- 
cording to Benetti et al. (2003), F(Ha) zvpc and A c (Ha) 
are also obtained from the radial ionization structure 
relative to ++ (Sect. 7.2) and the fair assumption 
0/H=constant across the nebula; at each position: 



F(Ha) z 



F(A5007A) 



zvpc 



cx | x /(T e ) x icf(0++), 



(9) 



with icf(0 ++ ) = opsonization correcting factor. 

For the internal, high-excitation regions icf(0 ++ ) 
comes from the ionization structure of helium (Seaton 
1968, Benetti et al. 2003): 



icf(0 ++ ) 



1 + 



0.3 x He ++ 



for the external, low-excitation layers we adopt: 

o° o+ 

icf(0++) outcr = l + -— + —-, (11) 



o++ 0++' 

which includes the ionization effects produced by the res- 
onant charge- exchange reaction + + H°=>0° + H + ; ac- 
cording to Osterbrock (1989) and Standi et al. (1999), in 
the outermost nebular regions H°/H+~0.83x (O /O+). 
Moreover (Sabbadin et al. 2004), 



N e (Ha) cx 



-0.47 



F(Ha) zvpc 1/2 
1 t| xD J ' 



(12) 



where D is the (known) distance, and ei the (un-known) 
"local filling factor" , i. e. the fraction of the local volume 
actually filled by matter with density N c . 

A comparative analysis gives quite satisfactory results: 
the iV e (Ha) profiles obtained in the two ways differ by 
less than 5% everywhere, except in the faint, innermost 
regions, where the discrepancy raises up to 10%. In the 
following we will adopt V c (Ha) given by F(A5007 A) zvpc 
and Eqs. (9) to (12). 

Fig. 13 shows the resulting radial distribution of the 
physical conditions in NGC 6741 at PA=15° and PA=75° 
(close to the apparent minor and major axes, respec- 
tively). 

The radial trend of T o [0 III] is common at all PA: 
>15 000 K in the weak innermost regions, rapidly decreas- 
ing outward down to 12 000 K in the densest layers, and 
more or less constant furtherout (uncertain). T C [N II] is 
systematically lower than T o [0 III]. All this is in quanti- 
tative agreement with the results by Hyung & Aller (1997) 
and Pottasch et al. (2001). 

N c presents a broad, asymmetric bell-shape profile 
with peaks up to 12 000 (±1000) cm~ 3 close to the mi- 
nor axis, and 6000-8000 (±800) cm -3 close to the major 
axis. Note that, in spite of the high density peaks - causing 
strong collisional de-excitation of the [S II] A6717 A line 
- we have N C [S II] ~ Ac (Ha), i. e. the local filling factor is 



ei~l, being A C [S II] 



xe 



0.5 



V (Ha) (Aller 1984, Pottasch 



He+ 



(10) 



1984, Osterbrock 1989). 

iV e [Ar IV] is systematically lower than iV c [S II] adopt- 
ing the electron impact excitation rates by Keenan et al. 
(1997), whereas N c [Ar IV] ~ N C [S II] for earlier collisional 
rates (for example, Aller 1984). 

Concerning the halo, N c peaks at the inner edge 
(~1000 - 1500 cm -3 ), and decreases outwards. 

Recent density determinations (mean values integrated 
over the slit) reported in the literature for NGC 6741 
are: ^(different ions)=6300 cm~ 3 by Hyung & Aller 
(1997), N C [S II]=6000 cm" 3 , N c [0 II]=8500 cm" 3 , N C [C\ 
III] =9000 cm~ 3 and N c [At IV] =5500 cm~ 3 by Pottasch 
ct al. (2001), and iV e [Cl III]=4470 cm" 3 , N c [At IV]=6610 
cur 3 and iV e [52/zrn/88/xm]=2880 cm" 3 by Liu et al. 
(2001). 

The ionized mass of NGC 6741, obtainable from the 
observed N c spatial distribution, the H/3 flux, and the ra- 
dio flux (Aller 1984, Pottasch 1984, Osterbrock 1989, and 
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Fig. 13. Radial distribution of the physical conditions (T e and N c ) at two selected directions of NGC 6741. Upper 
panel: the zvpc at PA=15° (close to the apparent minor axis); lower panel: the zvpc at PA=75° (close to the apparent 
major axis). Left ordinate scale: T c [0 III] (continuous line) and T [N II] (empty squares). Right ordinate scale: JV c (Ha) 
for D=2000 pc (dotted line), N C [S II] (empty circles), and iV e [Ar IV] (filled circles). 



Turatto et al. 2002), results to be M ion ~0.06(±0.02) M©, 
i. e. much lower than the mass of the external, neutral 
cocoon, M neutra i— 0.20 (±0.05) M Q (see Sect. 4), thus con- 
firming the peculiar evolutionary phase of our nebula. 

To be noticed: the relatively high N c of the halo implies 
a quite recent recombination start, no more that 200 years 
ago. In fact, the N c depletion rate for recombination is: 



diVc/dt = -as xN c x N(R+), (13) 



with aB=cffective recombination coefficient. Assuming 
N c ~ N(B+) and T c =12 000 K, and neglecting the recom- 
bination delay due to expansion, we obtain: 

N c (0) = JV e (t)/[l - 8.2 x 10- 6 x t x JV e (t)], (14) 

where N e (0) is the initial electron density, and N c (t) the 
electron density at time t (in years) elapsed from the re- 
combination start. Adopting 7V C (0)~10 000 cm~ 3 (mean of 
the density peaks in Fig. 13), and iV e (t)~600 cm -3 (the 
lowest value in Fig. 13), we have t~190 years. A very sim- 
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ilar "recombination age" is inferred from the mean H I 
density of the circum-nebular matter (Sect. 3). 

All this, combined with the short nebular age (~1400 
years), fairly agrees with the theoretical evolutionary 
times of a 0.66-0.68 M post-AGB star. According to 
Blocker (1995 and references therein), in the interval 
tsw (end of the superwind ejection) to tsw+H00yr the 
hydrogen-burning star crosses horizontaly the H-R dia- 
gram (at log L^/L© ~4.00) becoming hotter and hotter 
(up to logT* ~5.20 K). At tsw+HOOyr the nuclear fu- 
elling becomes insufficient, and in a century the luminosity 
declines to log L*/Lq ~3.70 (whereas T* reaches its maxi- 
mum, logT* ~5.30 K). At t S w+1200yr the hydrogen-shell 
burning ceases at all, the stellar evolution drastically ac- 
celerates till t sw +1400yr (log L„/L Q ~2.80; logT* ~5.20 
K), and later slows down (for example, at tsw+3000yr we 
have log L*/L Q ~2.40 and logT* ~5.15 K). 

Summing up: NGC 6741 is close to (or has even passed) 
the recombination-reionization transition. 

7.2. Radial ionic and mean chemical abundances, and 
photo-ionization model 

The zvpc of the different emissions furnishes the detailed 
radial ionization structure of NGC 6741. Since the blurred 
Ha appearance lowers the accuracy of the F(Ha) zvpc dis- 
tribution, we adopt A5007 A as reference emission, thus 
inferring the radial ionization structure relative to ++ 
(for details, see Benetti et al. 2003). The ^+ profiles of 
NGC 6741 at PA=15° and PA=75° (close to the apparent 
minor and major axes, respectively), presented in Fig. 14, 
confirm that the nebula is optically thick at all directions. 

Within this scenario, the [O III] rays punching the 
[N II] skin (Fig. 2) identify a few radial directions along 
and close to the true major axis of the nebula with "low" 
N e peaks (7V o (peak)<4000 cm" 3 according to the photo- 
ionization model presented in the course of the section): 
in these directions (barely visible in our spectral images at 
PA=75°, 95° and 115°; see Figs. 4 and 10) recombination 
processes are less efficient, and the gas can be (partially) 
ionized by the UV stellar flux at larger distanced. 

When combined with the results stored in the previous 
sections, all this indicates that the overall, complex struc- 
ture of NGC 6741, fully driven by the fast evolving cen- 
tral star, can be divided into three distinct zones showing 
quite different physical characteristics: the internal, high- 
to-medium excitation nebula directly ionized by the stellar 
radiation ends in a sharp transition zone at medium-to- 
low excitation marking the present edge of the ionization 
front. It is surrounded by the gas no more reached by 
the fading UV flux (i.e. the halo), dominated by recom- 
bination processes. A similar onion-like radial structure 
(ionization-transition-recombination) is observed in NGC 
6565 (Turatto et al. 2002). 

Though the ionization and thermal structure of NGC 
6741 result to be out of equilibrium, the large N e of the 
ionized gas (peaks up to 12 000 cm" 3 ) assures a quick 



nebular reply to the variable stellar flux. In particular, the 
delay in the "transition zone" is short, the recombination 
time, t rcc =l/(o!B x N c ), being a few years for hydrogen 
(even shorter for higher ionization species). 

In conclusion: NGC 6741 is nearly in equilibrium. This 
allows us to neglect (in first approximation) the time de- 
pendance, and estimate the total chemical abundances of 
the ionized gas with the classical method valid for "static" 
nebulae. 

We follow a two-step procedure: 
I (broad approach): according to the critical analy- 
ses by Alexander & Balick (1997) and Perinotto et al. 
(1998), the mean ionic abundances are obtained from the 
line fluxes integrated over the spatial profile and the ex- 
pansion velocity field; later on, we correct for the unob- 
served ionic stages by means of the ionization correcting 
factors, derived both empirically (Barker 1983, 1986), and 
from interpolation of theoretical nebular models (Shields 
et al. 1981, Aller & Czyzak 1983, Aller 1984, Osterbrock 
1989). The resulting mean chemical abundances are listed 
in Table 3 (last but one row); the recent estimates reported 
in the literature are also given in the Table. 

II (refining): we apply the photo-ionization code 
CLOUDY (Ferland et al. 1998) to a model-nebula charac- 
terized by the same distance, gas distribution, mean chem- 
ical composition and exciting star parameters of NGC 
6741, and then combine the model-line profiles (convolved 
for a seeing+guiding of 0.60") with the line profiles ob- 
served in the true nebula. 

Let us focus on the zvpc of the S-SW sector at PA=15° 
(close to the apparent minor axis of NGC 6741). The in- 
put parameters of the model-nebula are given in Table 
4, whereas Fig. 15 shows the physical conditions in the 
model-nebula (to be compared with Fig. 13), and the ab- 
solute radial flux distribution of the main emissions in 
both the model-nebula and NGC 6741. 

At first glance, Fig. 15 provides a satisfactory model- 
nebula vs. true-nebula agreement for the internal, high- 
to-medium excitation regions and the transition zone. 
Furthcrout, the weak recombining halo appears in NGC 
6741, whereas emissions are absent at all in the static 
model-nebula. The same is observed at the other PA. 

When entering in more detail, some disturbing drifts 
arise in Fig. 15. Two minor discrepacies are: 

(a) uncertain flux distribution of [Ne III] in NGC 6741 
(ascribable to the edge location of A3968 A in the orig- 
inal NTT+EMMI frames), 

(b) comparable model- nebula vs. true-nebula flux-shift 
for [S II] and [S III], indicating a slight underestimate 
of sulphur abundance in Table 3 (last but one row) , and 
in Table 4 (model- nebula). The same occurs for [N II]. 
The improved chemical abundances of NGC 6741 are 
reported in the last row of Table 3. 

The two major problems in Fig. 15 concern: 

(1) the ionization structure of helium: He I A5876 A and 
He II A4686 A are weaker and stronger, respectively, 
in the model-nebula than in the true-nebula; 
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Fig. 14. The radial ionization structure of NGC 6741 in the zvpc at PA=15° (close to the apparent minor axis; 
upper panel), and PA=75° (close to the apparent major axis; lower panel). Same orientation as Fig. 13. 



(2) the peak emission and whole flux distribution of the 
highest excitation ions (He II, [Ar IV] and [Ar V]) tend 
to be systematically closer to the central star in the 
true-nebula than in the model-nebula. 

Ad hoc manipulation of the model-nebula input data 
lowers - or even cancels - the foregoing problems; for exam- 



ple, we can act on the ionizing UV distribution, given the 
strong T» dependance of both the ionic structure of helium 
and the radial profile of the highest-excitation emissions. 
In fact, the only decrease by 17% of the stellar tempera- 
ture (i.e. for T* ~150000 K, which is within the observa- 
tional error box given in Sect. 6) simultaneously relieves 
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Reference 


He 


C 


N 


O 


Ne 


S 


CI 


Ar 


Hyung & Abler 1997 (icf) 


11.01 


8.86 


8.38 


8.73 


8.12 


6.89 


5.43 


6.49 


Hyung & Abler 1997 (model B) 


11.04 


8.90 


8.15 


8.65 


8.00 


6.76 


5.34 


6.54 


Pottasch et al. 2001 


11.04 


8.81 


8.45 


8.82 


8.26 


7.04 




6.69 


Pottasch et al. 2002 


11.04 


8.56 


8.26 


8.65 


8.18 


6.90 


5.26 


6.51 


Perinotto et al. 2004b 


11.08 




8.31 


8.70 


8.16 


6.70 




6.54 


This paper (icf) 


11.04 




8.20 


8.66 


8.08 


6.78 




6.54 


This paper (CLOUDY) 


11.04 




8.28 


8.66 


8.08 


6.90 




6.54 



Table 4. Input parameters for the model nebula 
(CLOUDY) 



Radial density profile 


ionized gas = Fig. 13; cfr. Sect. 7.1 




neutral gas = cfr. Sects. 3 and 4 


Chemical abundances: 




C, CI, K, Ca, Mg, Si 


Hyung & Aller (1997) model B 


He, N, O, Ne, S, Ar 


Table 3 (last but one row) 


other elements 


PN (CLOUDY default) 


Dust 


PN (CLOUDY default) 


Local filling factor 


1.0 




blackbody distribution 


Exciting star 


T» =170 000 K 




log L*/L Q = 2.75 


Distance 


2.0 kpc 



discrepancies (1) and (2). They almost disappear with a 
further, modest increase (by 20%) of the matter density 
in the innermost nebular layers. 

Unfortunately, the large number of involved parame- 
ters (density profile, chemical abundances, dust, temper- 
ature and luminosity of the star) and assumptions (dis- 
tance, black-body distribution, seeing convolution), added 
to the nebular compactness and peculiar evolutionary 
phase, greatly complicate the affair, whose definitive solu- 
tion needs detailed observations at even higher spatial and 
spectral resolutions, and is beyond the aims of the paper 
(as well beyond the patience of the authors). 



8. The 3-D morpho-kinematical structure 

According to the reconstruction method introduced by 
Sabbadin et al. (1985, 1987), and deepen by Sabbadin et 
al. (2000a, b), we select A4686 A of He II, A5007 A of [O 
III] and A6584 A of [N II] as markers of the high, medium 
and low- ionization regions of NGC 6741, respectively. The 
spectral images are: 

(a) de-convolved for seeing, spectral resolution and ther- 
mal motions (also fine-structure is taken into account for 
the recombination line of He II) , 

(b) de-projected through the relation V^ xp (km s _1 )= 
13 xR" (see Sect. 3), 

(c) assembled by means of the 3-D rendering procedure 
described by Ragazzoni et al. (2001). 

Here we present a limited number of frames, corre- 
sponding to a partial rotation around the North-South 
axis centered on the exciting star (i.e. almost perpendic- 
ular to the major axis). The complete series of nebular 
movies is available: 

- in the electronic version of the paper (on-line data), 

- at http://web.pd.astro.it/sabbadin the WEB 
site dedicated to the 3-D structure of expanding nebulae, 

- at http:/ /www. edpsciences.org 

Figs. 16 and 17 show the opaque reconstruction of 
NGC 6741 in He II, [O III] and [N II] at high and low flux- 
cuts, respectively, for a rotation of 180° through the first 
Euler angle. The (0,0) images correspond to the Earth- 
nebula direction (North is up and East to the left). An 
indication of the N c limit in Figs. 16 and 17 can be 
obtained from the corresponding [O III] absolute flux- 
cut: log E(A5007 A)= -15.95 and -16.70 erg s" 1 cm" 3 
(high and low cuts, respectively). Since O/H=4.6xl0 -4 
(Table 3), icf(0++)~l.l, and N c ~1.15xN(H+), they cor- 
respond to iV c (high cut)~9000 cm~ 3 for Fig. 16, and 
iV e (low cut)~4000 cm~ 3 for Fig. 17 (assuming T c =12 000 
K and ei=l). 
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MODEL NEBULA - PHYSICAL CONDITIONS 
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Fig. 15. Model-nebula (CLOUDY) vs. true-nebula (NGC 6741, PA=15°, S-SW sector): physical conditions and radial 
ionization structure. Top panel: physical conditions in the model-nebula (to be compared with Fig. 13). Second panel: 
absolute flux distribution of H I A6563 A, [N II] A6584 A, and the ionic sequence of helium (He I A5876 A and 
He II A4686 A) in the model- nebula (thin symbols) and the true-nebula (thick symbols). Third panel: absolute flux 
distribution for the ionic sequences of oxygen ([0 I] A6300 A, [O II] A7320 A, and [O III] A5007 A) and sulphur ([S II] 
A6717+6731 A, and [S III] A6312 A); same symbols as in the second panel. Bottom panel: absolute flux distribution 
for [Nc III] A3968 A, and the ionic sequence of argon ([Ar III] A7135 A, [Ar IV] A4740 A, and [Ar V] A7005 A); same 
symbols as in the second and third panels. 
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Fig. 16. Stereoscopic structure of NGC 6741 for a rotation around the North-South axis centered on the exciting 
star. Opaque reconstruction at high flux-cut for A4686 A of He II, A5007 A of [O III], and A6584 A of [N II], as seen 
from thirteen directions separated by 15°. The line of view is given by (0,tp), where 9 is the zenith angle and -0 the 
azimuthal angle. Each horizontal pair represents a "direct" stereoscopic pair (i. e. parallel eyes), and the whole figure 
provides twelve 3-D views of the nebula in as many directions, covering a straight angle (for details, see Ragazzoni et 
al. 2001). The (0,0) images represent the rebuilt-nebula as seen from Earth (North is up and East to the left). The 
complete scries of nebular movies is shown in the electronic version of the paper (on-line data). 
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The projection of NGC 6741 for a rotation around 
the N-S direction (almost perpendicular to the major 
axis) is presented in Fig. 18 (multi-color images only 
in the "free" electronic version of the paper, for lack 
of funds), providing a representative sample of the neb- 
ular appearance when changing the line of view. The 
left panel, (0,0), corresponds to NGC 6741 as seen from 
Earth (North is up and East to the left), to be compared 
with Fig. 1 and the HST image by Hajian & Terzian at 
http: / / ad. usno.navy.mil/pne/gallery.html 

The overall characteristics of NGC 6741 (a hot and 
faint star at the centre of a dense, inhomogeneous, equa- 
torial torus merging into a closed, ellipsoidal structure 
embedded into a large cocoon of almost neutral matter) 
closely resemble NGC 6565, a compact, recombining PN 
projected pole-on (Turatto el al. 2002), and show remark- 
able affinities with NGC 7027, the PNe prototype (Cox 
et al. 2002, Bains et al. 2003). Very likely, all three ob- 
jects represent "middle-aged" snapshots of a massive, fast 
evolving post-AGB star. 

9. Discussion 

In spite of the nebula compactness, we have quite success- 
fully applied the tomographic and 3-D analyses to NGC 
6741, covered at high spatial and spectral resolutions with 
ESO NTT+EMMI. 

NGC 6741 is a young (age~ 1400 years) PN at a dis- 
tance of 2000 pc; it consists of a high-excitation, almost- 
prolate ellipsoid (0.039 pc x 0.022 pc x 0.019 pc; major, 
intermediate and minor semi-axes, respectively) denser at 
the equator than at the poles (N Q peaks up to 12 000 
and 7000 cm -3 , respectively), surrounded by a thin skin 
at low-excitation (the transition zone, corresponding to 
the present ionization front), and embedded into an ex- 
tended, spherical (radius~0.080 pc), almost neutral re- 
combining halo containing a large fraction of the nebular 
mass (Mhaio ^0.20 M Q versus M ion ~0.06 M ). 

The complex ionization structure of NGC 6741 is fully 
driven by the fast evolution of the central star, a mas- 
sive (0.66-0.68 M Q ), hot (log T* ~5.23 K) and faint (log 
L*/Lq ~2.75) post-AGB star fading towards the white- 
dwarf region (after exhaustion of the hydrogen shell nu- 
clear burning). 

The ionized gas of NGC 6741 follows the general ex- 
pansion law T4 xp (km s _1 ) = 13xR", except the innermost, 
highest-excitation layers, showing clear evidences of decel- 
eration. This peculiar kinematical behaviour can be quali- 
tatively ascribed to the luminosity drop of the central star 
(Schonberner et al. 1997, Steffen et al. 1998, Marigo et al. 
2001, Pcrinotto et al. 2004a): in this evolutionary phase 
the stellar mass-loss quickly decreases, and the falling hot- 
bubble's pressure no more balances the pressure of the 
ionized gas, so that the contact discontinuity and the ad- 
jacent, high-excitation nebular layers accelerate inwards. 

The argument is poorly known, and deserves more at- 
tention, the kinematical properties of the innermost layers 
being an excellent diagnostic of the star-nebula interac- 



tion. A search for more "decelerated" candidates, a high- 
resolution spectroscopic survey in ionic species at large 
IP, and the detailed comparison with the current hydro- 
dynamical simulations are highly desired. 

Note that the opposite situation (i. e. acceleration) 
is expected for the external, recombining layers of NGC 
6741, whose de-pressuring gives rise to Rayleigh- Taylor 
instability at the ionization edge. Though the knotty and 
filamentary [N II] nebular appearance (Fig. 1 , lower panel) 
qualitatively agrees with this scenario, its confirmation 
needs very- high spectral resolution (R> 100 000) echel- 
lograms detecting the increase of expansion velocity (and 
turbulence!) at the ionization edge. 

A further interesting point concerns the wide radial 
matter distribution of NGC 6741 (much wider than the N e 
profile given in Fig. 13, due to the presence of the dense, 
almost-neutral recombining halo). This indicates a mod- 
est contribution of wind interaction to the nebular shaping 
(it essentially supports the innermost layers, avoiding gas 
infall), whereas the lion's share is taken by ionization: the 
thermal pressure of the ionized gas accelerates the mat- 
ter outwards and decelerates inwards, creating the density 
and velocity distributions observed in NGC 6741 (more 
general comments are towards the end of the section). 

Though the precise radial profile of the matter remains 
unknown - a large fraction of NGC 6741 being neutral -, 
a more sign provides that the outer radius of the main 
nebula extends well beyond the present ionization edge: 
the absence of any extended, diffuse, faint and roundish 
AGB-halo emission in the Palomar Sky Survey plates, 
deep imaging by Schwarz et al. (1992), and HST frames. 
This is confirmed by our spectra: no kinematical signature 
of an AGB-halo (i. e. external, faint, un-tilted and un-split 
emission) appears. Since the low-density gas of an AGB- 
halo (usually Ac=50 to 200 cm -3 ; Chu et al. 1987; Corradi 
et al. 2003) is little affected by recombination (from Eq. 
(14), Ac(0)=50 to 200 cm- 3 implies AT e (200yr)=47 to 162 
cm -3 ), we infer that: 

- NGC 6741 never became optically thin to the UV stellar 

radiation, 

- the total nebular mass, M tot , is larger than 0.40 M Q . 

Let us focus on recombination, which appears as the main 
characteristic (and peculiarity) of our nebula. As outlined 
in the previous sections, NGC 6741 is the second PN - 
of the six so far analysed with the 3-D methodology - in 
a deep recombination phase (the first being NGC 6565, 
Turatto et al. 2002). Moreover, according to Benetti et al. 
(2003), a third nebula, NGC 6818, is at the very beginning 
of recombination. Such an apparent excess of recombining 
PNe is a mere consequence of our target selection crite- 
ria (in particular, high surface brightness and richness of 
ionic species), favouring massive nebulae powered by fast 
evolving central stars. 

Thanks to the unprecedented accuracy achieved by 
tomography and 3-D recovery on the manifold facets of 
the PN research, we have identified a number of observa- 
tional evidences supporting the recombination hypothesis 
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for NGC 6741; they are: 

(1) nebula much brighter than the star (Sects. 1 and 6), 

(2) presence of absorbing knots (Sect. 1), 

(3) co-existence of ionic species in a large IP range (Sects. 
1, 3 and 7.2), 

(4) kinematics of the halo (Sect. 3), 

(5) high density of the almost-neutral halo (Sects. 3 and 

4), 

(6) c(H/3) variable over the nebula (Sect. 4), 

(7) high temperature and low luminosity of the central 
star (Sect. 6), 

(8) T Z H I>T z Hc II (Sect. 6). 

A minimal, but effective, recombination standard for 
the whole PN class - in practice, a summary of points (1) 
and (7) - is given by: 

m(V,) obs + logF(H/3) obs - 1.1 x c(H/3) > 5.0, (15) 

where F(H/3) b s is in erg cm~ 2 s _1 . 

A quick look at the Strasbourg-ESO Catalogue of 
Galactic PNe (Acker et al. 1992) allowed us to select some 
fifty targets satisfying Eq. (15). Half of them (e. g. NGC 
2440, Hb 4, NGC 6565, NGC 6537, NGC 6620, NGC 6741, 
NGC 6886, NGC 6881, NGC 7027, Hu 1-2, NGC 6302, 
and NGC 6563) are compact (sometimes with butterfly- 
like cxtcntions), at high surface brightness, and represent 
young recombining PNe ejected and excited by a massive 
(M*>0.64 M ) fast evolving star beyong the turn-around 
point in the H R diagram. The remaining candidates are 
mean-to-low surface brightness, extended (often bipolar) 
nebulae (e. g. NGC 6445, NGC 6439, A 55, NGC 6772, A 
53, NGC 6818, NGC 6894, NGC 7048, NGC 2438, NGC 
2818, and IC 4406), corresponding to aged, evolved PNe, 
optically thick (at least in some directions) to the fading 
UV radiation of a moderately-massive (0.60 to 0.64 M©) 
post-AGB star in the white dwarf cooling sequence. 

Further supports to the "massive" nature for the cen- 
tral stars of PNe satisfying Eq. (15) are: 

- most targets belong to the morphological class B of 

Grcig (1971, 1972), 

- a large fraction consists of Type-I PNe (Peimbert 1978; 

Pcimbcrt & Torres-Peimbert 1983; Torres-Peimbert & 
Peimbert, 1983; Peimbert et al. 1995), 

- the H2 emission at 2.12 /im is commonly detected 

(Kastner et al. 1996; Natta & Hollenbach 1998; 
Bohigas 2001). Note that, though NGC 6741 is not 
listed among the H 2 -emitters, a deep search promises 
fruitful results. 

The general rule for recombining PNe is that, the 
higher the surface brightness of the nebula, the larger is 
the mass of the star. 

Summing up: recombination represents a not un- 
common evolutionary phase in PNe; it becomes un- 
escapable for massive nebulae ejected and excited by a 
massive, fast-evolving post-AGB star (also see Tylenda 
1986, Szczerba 1990, Phillips 2000, Marigo et al. 2001, 
Turatto et al. 2002, Benetti et al. 2003, and Perinotto et 
al. 2004a). 



Let us close this verbose paper with some general con- 
siderations on the kinematics of the whole PN class. The 
comparative analysis performed in Sect. 3 - rejecting the 
kinematical results based on spherical symmetry assump- 
tion and/or emission profiles of recombination lines inte- 
grated along the slit - enhances the rarity of real nebulae 
showing a "U'-shaped expansion velocity field. 

We stress that: 

(a) Wilson's law (V exp oc radius) is the general rule for 
PNe (Wilson 1950, Weedman 1968, Sabbadin et al. 
2004 and references therein); 

(b) at present, only BD+30°3639 and NGC 40 show clear 
observational evidences of a soft acceleration in the 
innermost layers (Bryce & Mellema 1999; Sabbadin 
et al. 2000a). Both objects: (I) are very-low excita- 
tion PNe (I([0 III] A5007 A)< I(H/3)), (II) exhibit an 
extremely sharp radial matter profile (shell thickness 
AR/RX0.25), and (III) are powered by a "cold" and 
bright central star of late- WR spectral type. According 
to Iben et al. (1983), Blocker (1995, 2001), and Herwig 
et al. (1999), they represent "born-again" PNe (i. e. 
the hydrogen-deficient star suffers a late thermal pulse 
during the motion towards the white-dwarf region. Due 
to the pulse-induced convection, hydrogen mixes and 
burns on a convective turn-over time scale); 

(c) three more "born-again" candidates - A 30, A 58 and 
A 78 - are faint, extended PNe presenting a central, 
complex structure of fast, hydrogen-deficient knots 
(Jacoby 1979, Pollacco et al. 1992). 

All this calls for severe constraints to the evolutive pa- 
rameters responsible of PNe shape and shaping (in par- 
ticular, wind interaction). In fact, according to the de- 
tailed radiation-hydrodynamics simulations (Perinotto et 
al. 2004a, and references therein), a PN is the result of 
the synergistic effects of ionization and fast wind on the 
gas ejected during the superwind phase, generating a typ- 
ical double-shell structure (inner "rim" + outer "shell") 
characterized by a "U" or better "V" -shaped expansion 
velocity field. The double-shell structure may only be de- 
stroyed either by recombination of the "shell" when the 
central star fades, or by overtaking of the "shell" by the 
"rim" . In both cases a single-shell configuration emerges, 
whose expansion velocity field is simple and increases al- 
ways steadily inwards. 

Wind interaction being the main responsible of the 
"V" -shaped expansion profile in double-shell model-PNe, 
as well of the increasing inwards velocity field in single- 
shell model-PNe, we infer that all current radiation- 
hydrodynamics simulations tend to overestimate 
the dynamical effects of fast stellar wind on the 
nebular gas (we suspect that post-AGB mass-loss rates 
adopted by the theoretical simulations are systematically 
too high). Though the same suspect arose in the study of 
NGC 6565 (Turatto et al. 2002), NGC 6818 (Benetti et 
al. 2003) and NGC 7009 (Sabbadin et al. 2004), a quanti- 
tative answer (also including: (a) the temporal evolution 
of mass-loss in the superwind phase, Schonberner et al. 
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2005, and (b) the binarity and the possible role of mag- 
netic fields, Garcia-Segura & Lopez 2000, Blackman et al. 
2001) will come from the detailed analysis performed on 
a representative sample of PNe in both hemispheres, cov- 
ered with ESO NTT+EMMI and TNG+SARG (work in 
preparation) . 

10. Conclusions 

In our opinion, PN observers should cancel the word "av- 
erage" (and respective synonyms) from their scientific vo- 
cabulary: no more average line fluxes, integrated spectral 
profiles, mean electron temperature, overall electron den- 
sity, and so on. In fact, to bridge the historical gap be- 
tween theory and practice (Aller 1994), we need detailed, 
point-to-point flux and velocity measurements in a wide 
ionization range, combined with a straightforward, versa- 
tile methodology of analysis. To this end, we have applied 
tomography and 3-D recovery (Sabbadin et al. 2004 and 
references therein) to long-slit echellograms of the com- 
pact, bright PN NGC 6741, covered at nine PA with ESO 
NTT+EMMI. 

Our long-lasting journey started with the gas kine- 
matics, penetrated into the galactic and circum-nebular 
absorptions, and explored the nebular distance, mass and 
age. Later on, we dealt with the stellar properties, pushed 
through the nebular physical conditions, ionization struc- 
ture and photo-ionization model, went over the stumbling 
block of image de-projection, and, at last, flew about the 
multi-color and opaque reconstructions of the nebula. 

The wealth and accuracy of the results here presented: 

- confirm the unique capacity of tomography and 3-D re- 

covery in extracting the huge amount of physical infor- 
mation stored in high-resolution spectroscopy (a tacit 
invitation, extended to the worldwide researchers), 

- stress once more that properly squeezed echellograms 

do represent a peerless tool for deeping the kinemat- 
ics, physical conditions, ionization structure and evolu- 
tion of all classes of expanding nebulae (PNe, nova and 
supernova remnants, shells around Population I Wolf- 
Rayet stars, nebulae ejected by symbiotic stars, bub- 
bles surrounding early spectral-type Main Sequence 
stars etc.). 
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Fig. 17. Same as Fig. 16, but at low flux-cut. 
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Fig. 18. Optical appearance of NGC 6741 (multi-color in the electronic version of the paper; bluc=Hc II, green=[0 III], 
red=[N II]) for a rotation through the N-S axis centered on the exciting star. The left panel, (0,0), corresponds to the 
re-built nebula as seen from the Earth (North is up and East to the left). Recall that projection^, ^)=projection(# ± 
180°, ip ± 180°). The complete (multi-color) movie is shown in the electronic version (on-line data). 



